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Abstract 



The resistance is computed of an NI1NI2S junction, where N = normal metal, 
S = superconductor, and Ij = insulator or tunnel barrier (transmission prob- 
ability per mode Tj). The ballistic case is considered, as well as the case 
that the region between the two barriers contains disorder (mean free path /, 
barrier separation L) . It is found that the resistance at fixed r2 shows a min- 
imum as a function of Ti, when Fi ~ ^/2T2, provided / ^ T2L. The minimum 
is explained in terms of the appearance of transmission eigenvalues close to 
one, analogous to the "reflectionless tunneling" through a NIS junction with a 
disordered normal region. The theory is supported by numerical simulations. 



I. INTRODUCTION 

Reflectionless tunneling is a novel quantum interference effect which occurs when dissi- 
pative normal current is converted into dissipationless supercurrent at the interface between 
a normal metal (N) and a superconductor (S) |T|. Experimentally, the effect is observed as 
a peak in the differential conductance around zero voltage or around zero magnetic field . 
Its name refers to the fact that, for full phase coherence, the Andreev-refiected quasiparticle 
can tunnel through the potential barrier at the NS interface without suffering reflections. 
(The potential barrier can be the insulator (I) in an NIS junction, or the Schottky barrier 
in a semiconductor-superconductor junction.) Application of a voltage or magnetic field 
destroys the phase coherence between electrons and holes, and thus reduces the conduc- 
tance of the junction. We now have a good theoretical understanding of the effect, based 
on a combination of numerical and analytical work 0-0. The basic requirement for 
reflectionless tunneling is that the normal region has a resistance which is larger than the 
resistance of the interface. In that case the disorder is able to open a fraction of the tunnel- 
ing channels, i.e. it induces the appearance of transmission eigenvalues close to one [0. As 
a result of these open channels, the resistance has a linear dependence on the transparency 



of the interface, instead of the quadratic dependence expected for Andreev reflection fTl 
(which is a two-particle process). 

The purpose of this work is to present a study of reflectionless tunneling in its simplest 
form, when the resistance of the normal metal is due to a second tunnel barrier, in series 
with the barrier at the NS interface. This allows an exact calculation, which shows many of 
the features of the more complicated case when the resistance of the normal region is due to 
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disorder. Furthermore, the double-barrier geometry provides an experimentally realizable 
model system, for example in tunneling from an STM into a superconductor via a metal 
particle . 

The outline of this paper is as follows. In Section 2 we consider the problem of a NI1NI2S 
junction without disorder. We compute the resistance of the junction as a function of the 
transmission probabilities per mode Fi and F2 of the two barriers. The resistance at fixed 
F2 shows a minimum as a function of Fi when Fi ^ V^r2 = T- The resistance in the 
minimum depends linearly on 1/F, in contrast to the quadratic dependence in the case of 
a single barrier. In Section 3 we apply a recent scaling theory [^, to find the influence 
on the resistance minimum of disorder in the region between the barriers (length L, mean 
free path /). The resistance minimum persists as long as / ^ TL. In the diffusive regime 
(/ ^ L) our results agree with a previous Green's function calculation by Volkov, Zaitsev, 
and Klapwijk 0. The analytical results are supported by numerical simulations, using the 
recursive Green's function technique ||13| . We conclude in Section 4. 



II. NINIS JUNCTION WITHOUT DISORDER 

We consider a NI1NI2S junction, where N = normal metal, S = superconductor, and Ij 
= insulator or tunnel barrier (see inset of Fig. 1). The transmission probability per mode of 
Ij is denoted by F,. For simplicity we neglect the mode- dependence of Fj. In this section, we 
assume ballistic motion between the barriers. (The effect of disorder in the normal region is 
considered in Sec. III.) A straightforward calculation yields the transmission probabilities 
Tn of the two barriers in series, 

Tn = {a + b cos (j)n)~^ , (2.1) 

where 

2 - Fi - F2 , , 

a = 1 + (2.2a) 



riF2 

2(l-F0V^(l-F2)^/^ 

riF2 



(2.2b) 



and (pn is the phase accumulated between the barriers by mode n = 1,2, ■ ■ ■ N (with the 
number of propagating modes at the Fermi level). If we substitute Fj = l/cosh^Oj (oj > 0), 
the coefficients a and b can be rewritten as 

a = I + I cosh 2a;i cosh 2q;2, (2.3a) 
6=1 sinh 2ai sinh 2a2. (2.3b) 

Since the transmission matrix t is diagonal, the transmission probabilities T„ are identical 
to the eigenvalues of tt^. We use the general relationship between the conductance Gns = 
Gninis of the NINIS junction and the transmission eigenvalues of the normal region 
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p{T) = ^-{h'T^-{aT-lfy"\ (2.8) 



which is the analogue of the Landauer formula, 

G'n = ^ E Tn, (2.5) 

for the conductance Gn = Cninin in the normal state. We assume that L ^ Af (Af is 
the Fermi wavelength) and NTi ^ 1, so that the conductance is not dominated by a single 
resonance. In this case, the phases 0„ are distributed uniformly in the interval (0, 2tt) and 
we may replace the summations in Eqs. ( |2.4| ), ( p.5| ) by integrals over 0: Y.n=i fi'Pn) 
{N / 2%) J^^ d(pf {(/)). The resuh is 

Ae^N cosh 2ai cosh 2^2 

Gns = — 7 — -3725 (2-6) 

^ (cosh^ 2ai + cosh^ 2^2 - l) 

4:6^ N 

Gn = — (cosh2ai + cosh 2^2)" • (2.7) 
h 

These expressions are symmetric in the indices 1 and 2: it does not matter which of the two 
barriers is closest to the superconductor. 

In the same way we can compute the entire distribution of the transmission eigenvalues, 
P{T) = EnS{T - Tn) ^ (Ar/27r)/(f"d05(T - T(0)). Substituting T(0) = (a + fccos*^)"^ 
from Eq. ( p.!] ), we find 

vrT 

In Fig. 1 we plot the resistance = 1/Gn and i?NS = V^nS) following from Eqs. ( p.6|) 
and (]2.7|). Notice that follows Ohm's law, 

R^ = ^(y^^ + '^/^^''^)^ (2-9) 

as expected from classical considerations. In contrast, the resistance i?Ns has a minimum 
if one of the F's is varied while keeping the other fixed. This resistance minimum cannot 
be explained by classical series addition of barrier resistances. If F2 -C 1 is fixed and Fi is 
varied, as in Fig. 1, the minimum occurs when Fi = v^F2. The minimal resistance R^s^ is 
of the same order of magnitude as the resistance i?N in the normal state at the same value of 
Fi and F2. (For F2 <^ 1, -Rn" = 1.52 _Rn) In particular, we find that -Rn™ depends linearly 
on 1/Fj, whereas for a single barrier i?NS oc 1/F^. 

The linear dependence on the barrier transparency shows the qualitative similarity of 
a ballistic NINIS junction to a disordered NIS junction. To illustrate the similarity, we 
compare in Fig. 2 the densities of transmission eigenvalues throught the normal region. The 
left panel is for an NIS junction (computed using the results of Ref. P]), the right panel is 
for an NINIS junction (computed from Eq. (|2.8|) ). In the NIS junction, disorder leads to a 
bimodal distribution p(T), with a peak near zero transmission and another peak near unit 
transmisssion (dashed curve). A similar bimodal distribution appears in the ballistic NINIS 
junction, for approximately equal transmission probabilities of the two barriers. There are 
also differences between the two cases: The NIS junction has a unimodal p{T) if L/l < 1/F, 
while the NINIS junction has a bimodal p{T) for any ratio of Fi and F2. In both cases, the 
opening of tunneling channels, i.e., the appearance of a peak in p(T) near T = 1, is the 
origin for the 1 /F dependence of the resistance. 
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III. EFFECTS OF DISORDER 



Let us now investigate what happens to the resistance minimum if the region of length L 
between the tunnel barriers contains impurities, with elastic mean free path /. We denote s = 
L/l. When introducing disorder, it is necessary to consider ensemble-averaged quantities. 
To calculate the ensemble-averaged conductance (Gns); we need to know the density p 
of the transmission eigenvalues T„ as a function of s. It is convenient to work with the 
parameterization 

Tn = 1/ cosh^Xn, x„ > 0. (3.1) 

The density of the x„'s is defined by p(x, s) = {J2n ^{x~Xn))- From Eq. ( ^.1| ) we know that, 
for s = (no disorder). 



f^'^ d(f) I 

p(x, 0) = N —— S{x — arccoshA/ a + b cos < 



— sinh 2x (b'^ - (a - cosh^ xf) ' , (3.2) 



TT 



for arccoshi/ct — b = a^min < a; < Xmax = arccosh-^/ a + b. 

For s > we obtain the density p{x, s) from the integro-differential equation [15 



-^p(x,s) = 'K~p{^^s)'K~ f dx'p(x', s) In I sinh^ X — sinh^ x'l, (3.3) 

OS 2N ox ox Jo 



which is the large A^-limit of the scaling equation due to Dorokhov |16] and Mello, Pereyra, 



and Kumar ||T^. This equation describes the evolution of p{x, s) when an infinitesimal slice 
of disordered material is added. With initial condition ( p.2|) it therefore describes a geometry 
where all disorder is on one side of the two tunnel barriers, rather than in between. In fact, 
only the total length L of the disordered region matters, and not the location relative to 



the barriers. The argument is similar to that in Ref. |]18[. The total transfer matrix M 
of the normal region is a product of the transfer matrices of its constituents (barriers and 
disordered segments): M = M1M2M3 ■ ■ -. The probability distribution of M is given by the 
convolution p{M) = pi * p2 * * ■ ■ ■ of the distributions Pi of transfer matrices Mj. The 
convolution is defined as 

p-*p.{M)= J dM'pi{MM'-^)pj{M'). (3.4) 

If for all parts i of the system, Pi{Mi) is a function of the eigenvalues of M^mJ only, the 
convolution of the pi commutes |l^. The distributions pi are then called isotropic. A 
disordered segment (length L, width W) has an isotropic distribution if L 3> W. A planar 
tunnel barrier does not mix the modes, so a priori it does not have an isotropic distribution. 
However, if the mode- dependence of the transmission probabilities is neglected (as we do 
here), it does not make a difference if we replace its distribution by an isotropic one. The 
commutativity of the convolution of isotropic distributions implies that the location of the 
tunnel barriers with respect to the disordered region does not affect p{x, s). The systems in 
Figs. 3a, b, and c then have identical statistical properties. 
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Once p(x, s) is known, the conductances (Gns) and (Gn) can be determined from 



(G'ns) 
(Gn) 



4g2 ^oo S 

— — / dx 2~ 

2e^ /"°°ci 
/i Jo cosh^ X 



cosh 2a; 
p{x,s) 



(3.5) 
(3.6) 



where we have substituted Eq. ( p.l|) into Eqs. (|2.4| ), (|2.5|) . In Ref. p a general solution 
to the evolution equation was obtained for arbitrary initial condition. It was shown that 
Eq. ( p.3|) can be mapped onto Euler's equation of hydrodynamics 



-U{C,s) + U{C,s)—U{C,s) 



0, 



(3.7) 



by means of the substitution 



^(C,^) 



sinh2C f°° 



2N 



p{x',s) 



sinh^ ( — sinh^ x' 



(3.8) 



Here, U = Ux + iUy and ( = x + iy. Eq. ( p.7| ) describes the velocity field U{(, s) of a 2D 
ideal fluid at constant pressure in the x-y plane. Its solution i^ 



U{C,s) = Uo{C-sU{C,s)), 



(3.9) 



in terms of the initial value Uq{Q = f/(C,0). The probability distribution p{x,s) follows 
from the velocity field by inversion of Eq. (|3.8|), 



p{x,s) 

where e is a positive infinitesimal. 



2N 



vr 



Uy{x - ie, s), 



(3.10) 



In our case, the initial velocity field [from Eqs. (3^) and 
Uo{C) = -\ sinh2C [(cosh^ C - af - 



IS 



(3.11) 



The resulting density (|3.1CI|) is plotted in Fig. 4 for Fi = r2 = F and several disorder 
strengths. The region near x = is of importance for the conductance (since x near zero 
corresponds to near-unit transmission). The number A'^open = p(0, s) is an estimate for 
the number of transmission eigenvalues close to 1 (so-called "open channels" |ll9[). In the 
absence of disorder, iVopen is non-zero only if Fi ^ F2 (then a — h = 1 ^ Xmin = 0). From 
Eq. (|3.2|) we find Nopen = NT /it for s = and Fi = F2 = F ^ 1. Adding disorder reduces 
the number of open channels. If Fi 7^ F2 there are no open channels for s = (xmin > 0). 
Disorder then has the effect of increasing A^open, such that A'opon ~ N/s if (Fi + F2)s ^ 1. 



^ The implicit equation (|3.9| ) has multiple solutions in the entire complex plane; we need the 
solution for which both ( and C — sU{Ci s) lie in the strip between the lines y = and y = —7r/2. 
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The disorder-induced opening of channels was studied in Refs. [0,|TO| for the case of a single 
tunnel barrier. 

To test our analytical results for the eigenvalue density p{x,s), we have carried out 
numerical simulations, similar to those reported in Ref. 0. The sample was modeled by a 
tight-binding Hamiltonian on a square lattice with lattice constant a. The tunnel barriers 
were accounted for by assigning a non-random potential energy Ub = 2.3 Ep to a single row of 
sites at both ends of the lattice, which corresponds to a mode averaged barrier transparency 
Fi = r2 = 0.18. The Fermi energy was chosen at 1.5 Mq, with uq = fi^/2ma?'. Disorder 
was introduced by randomly assigning a value between ±|^7d to the on-site potential of 
the lattice points between the barriers. The disorder strength Ud was varied between 
and 1.5 Mo, corresponding to s between and 11.7. We considered geometries with both 
a square disordered region (285 x 285 sites, = 119) and a rectangular one (285 x 75 
sites, N = 31), to test the geometry dependence of our results. In Fig. 5, we compare the 
integrated eigenvalue density z/(x, s) = Jq dx'p{x', s) with the numerical results. 

The quantity z/(a;, s) follows directly from our simulations, by plotting the x„'s in as- 
cending order versus n/N = v. We want to sample v{x, s) at many points along the x-axis, 
so we need N large. Since the self-averaging (fluctuations are of the order of 1/iV), 

it is not necessary to average over many samples. The data shown in Fig. 5 are from a 
single realization of the impurity potential. There is good agreement with the analytical 
results. No geometry dependence is observed, which indicates that the restriction L ^ W 
of Eq. (|3.3| ) can be relaxed to a considerable extent. 

Using Eqs. (|3.5| ) and (|3.^ ), the average conductance (Gns) can be directly expressed in 
terms of the velocity field, 

2Ne^ d 

(G,,)_Jta^_f/(C..). (3.12) 
For C — ivr/4, U — * iUy, Uy > 0. The implicit solution (|3.9|) then takes the form 



0^(2a + sin0- 1)2-462 = 2scos(j), (3.13) 

where = 2sUy G [0,7r/2]. We now use that ^U{C,s)\^^_^ = -[£f/(-f , s)]/f/(-if , s) 
[see Eq. (|3)]. Combining Eqs. (|3l2|) and ( gll we find 

9A^(=2 

(Gns) = -^(3 + 1/Q)-\ (3.14) 
where the effective tunnel rate Q is given in terms of the angle (f) in Eq. ( p.l3| ) by 

Q = 77^ fsin + ^ [sin + (1 - 2/ri)(l - 2/r2)]^ . (3.15) 



SCOSi 



Eqs. ( ^.13| )- p.l5| ) completely determine the conductance of a double-barrier NS-junction 
containing disorder. 

In Fig. 6, we plot (-Rns) for several values of the disorder, keeping r2 = 0.1 fixed and 
varying the transparency of barrier 1. For weak disorder (r2S <^ 1), the resistance minimum 
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is retained, but its location moves to larger values of Fi. On increasing the disorder, the 
minimum becomes shallower and eventually disappears. In the regime of strong disorder 
(r2S ^ 1), the resistance behaves nearly Ohmic. 

We stress that these results hold for arbitrary s = L/l, all the way from the ballistic into 
the diffusive regime. Volkov, Zaitsev, and Klapwijk ||^ have computed (Gns) in the diffusive 
limit s ^ 1. In that limit our Eqs. ( p. 131) and (|3.15|) take the form: 



SCOS( 



1 



1 



Tos 



1 



1 



(3.16) 



1 

Q 



-1/2 



TiS 



(3.17) 



in precise agreement with Ref. 0. Nazarov's circuit theory ||T0[, which is equivalent to the 
Green's function theory of Ref. 0, also leads to this result for (Gns) in the diffusive regime. 

Two limiting cases of Eqs. ( |3.16| ) and ( p.l7| ) are of particular interest. For strong barriers, 
Fi, r2 -C 1, and strong disorder, s 3> 1, one has the two asymptotic formulas 



(Gns) 
{Gns) 



p2p2 



h (rf + r 



2\3/2' 
2) 



if ri,r2<i/s. 



h 



S + l/Fi + 1/12 



if ri,r2>i/s. 



(3.18) 
(3.19) 



Eq. ( p. 181 ) coincides with Eq. ( ^.61 ) in the limit ai,a2 ^ 1 (recall that Fj = 1/ cosh 



a,; 



This shows that the effect of disorder on the resistance minimum can be neglected as long 
as the resistance of the junction is dominated by the barriers. In this case (Gns) depends 
linearly on Fi and F2 only if Fi ^ F2. Eq. ( |3.19| ) shows that if the disorder dominates, (Gns) 
has a linear F-dependence regardless of the relative magnitude of Fi and F2. 



IV. CONCLUSIONS 

In summary, we have derived an expression for the conductance of a ballistic NINIS 
junction in the limit A^F ^ 1 that the tunnel resistance is much smaller than h/e^. In this 
regime the double-barrier junction contains a large number of overlapping resonances, so that 
in the normal state the resistance depends monotonically on 1/F. In contrast, the NINIS 
junction shows a resistance minimum when one of the barrier transparencies is varied while 
the other is kept fixed. The minimal resistance (at Fi ~ F2 = F) is proportional to 1/F, 
instead of the 1/F^ dependence expected for two-particle tunneling into a superconductor. 
This is similar to the reflectionless tunneling which occurs in an NIS junction. Using the 
results of the ballistic junction, we have described the transition to a disordered NINIS 
junction by means of an evolution equation for the transmission eigenvalue density 0. We 
found that the resistance minimum is unaffected by disorder, as long as I ^ L/T, i.e., as long 
as the barrier resistance dominates the junction resistance. As the disorder becomes more 
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dominant, a transition to a monotonia F- dependence takes place. In the limit of diffusive 
motion between the barriers, our results agree with Ref. 0. 

Throughout this paper we have assumed zero temperature, zero magnetic field, and 
infinitesimal applied voltage. Each of these quantities is capable of destroying the phase 
coherence between the electrons and the Andreev-reflected holes, which is responsible for 
the resistance minimum. As far as the temperature T and voltage V are concerned, we 
require kBT,eV -C h/rdweii for the appearance of a resistance minimum, where Xdweii is the 
dwell time of an electron in the region between the two barriers. For a ballistic NINIS 
junction we have Xdweii ~ L/vpT, while for a disordered junction Td^eii ~ L'^/vpTl is larger 
by a factor L/l. It follows that the condition on temperature and voltage becomes more 
restrictive if the disorder increases, even if the resistance remains dominated by the barriers. 
As far as the magnetic field B is concerned, we require B <^ h/eS (with S the area of the 
junction perpendicular to B), if the motion between the barriers is diffusive. For ballistic 
motion the trajectories enclose no flux, so no magetic field dependence is expected. 

A possible experiment to verify our results might be scanning tunneling microscopy of 
a metal particle on top of a superconducting substrate |T2|. The metal-superconductor 
interface has a fixed tunnel probability F2. The probability Fi for an electron to tunnel 
from STM to particle can be controlled by varying the distance. (Volkov has recently 
analyzed this geometry in the regime that the motion from STM to particle is diffusive 
rather than by tunneling [^.) Another possibility is to create an NINIS junction using a 
two-dimensional electron gas in contact with a superconductor. The tunnel barriers could 
then be implemented by means of two gate electrodes. In this way both transparancies 
might be tuned independently. 
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Figure 1 

Dependence of the resistances i?N and i?NS of ballistic NININ and NINIS structures, respec- 
tively, on barrier transparancy Fi, while transparancy r2 = 0.1 is kept fixed [computed from 
Eqs. (|2.6| ) and (|2.7| )]. The inset shows the NINIS structure considered. 
Figure 2 

Density of normal-state transmission eigenvalues for an NS junction with a potential barrier 
at the interface (transmission probability F = 0.4). The left panel (a) shows the disorder- 
induced opening of tunneling channels in an NIS junction (solid curve: s = 0.04; dotted: 
s = 0.4; dashed: s = 5; where s = L/l). The right panel (b) shows the opening of channels 
by a second tunnel barrier (transparancy V) in an NINIS junction (solid curve: V = 0.95; 
dotted: T' = 0.8; dashed: T' = 0.4). The curves in (a) are computed from Ref. the 
curves in (b) from Eq. ( p.8| ). Notice the similarity of the dashed curves. 
Figure 3 

The systems a,b, and c are statistically equivalent, if the transfer matrices of each of the 
two barriers (solid vertical lines) and the disordered regions (shaded areas, Li + L2 = L in 
case b) have isotropic distributions; in that case, the position of the disorder with respect 
to the barriers does not affect the eigenvalue density p(x, s). 
Figure 4 

Eigenvalue density p(x, s) as a function of x (in units of s = L/l) for Fi = r2 = 0.2. Curves 
a,b,c,d, and e are for s = 0.5, 2, 5, 20, 100, respectively. In the special case of equal tunnel 
barriers, open channels exist already in the absence of disorder. 
Figure 5 

Comparison between theory and simulation of the integrated eigenvalue density for Fi = 
F2 = 0.18. The labels a,b,c indicate, respectively, s = 0,3,11.7. Sohd curves are from 
Eq. (^.91); data points are the x„'s from the simulation plotted in ascending order versus 
n/N. Filled data points are for a square geometry, open points are for an aspect ratio 
L/W = 3.8. 
Figure 6 

Dependence of the ensemble-averaged resistance (-Rns) for a disordered NINIS junction on 
barrier transparancy Fi, while F2 = 0.1 is kept fixed [computed from Eqs. (|3.14| ) and p.l5| )]. 
Curves a, b, c, d are for s = 0, 2, 7, 30, respectively. The resistance minimum persists for small 
disorder. 
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